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Given oracle access to an unknown unitary C from the Clifford group and its conjugate, we give 
an exact algorithm for identifying C with O(n) queries, which we prove is optimal. We then extend 
this to all levels of the Gottesman-Chuang hierarchy (also known as the Ck hierarchy). Further, for 
unitaries not in the hierarchy itself but known to be close to an element of the hierarchy, we give 
a method of finding this close element. We also present a Clifford testing algorithm that decides 
whether a given black-box unitary is close to a Clifford or far from every Clifford. 
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I. INTRODUCTION 

A central problem in quantum computing is to de- 
termine an unknown quantum state from measurements 
of multiple copies of the state. This process is known 
as quantum state tomography (see [1] and references 
therein). By making enough measurements, the prob- 
ability distributions of the outcomes can be estimated 
from which the state can be inferred. A related prob- 
lem is that of quantum process tomography, where an 
unknown quantum evolution is determined by applying 
it to certain known input states. There are several meth- 
ods for doing this, including what are known as Standard 
Quantum Process Tomography [2, 3] and Ancilla Assisted 
Process Tomography [4, 5]. These methods work by us- 
ing state tomography on the output states for certain 
input states. 

However, all these procedures share one important 
downside: the number of measurements required in- 
creases exponentially with the number of qubits. This 
already presents problems even with systems achievable 
with today's technology, for which complete tomographi- 
cal measurements can take hours (e.g. [6]) making tomog- 
raphy of larger systems unfeasible. Unfortunately this 
exponential cost is necessary to determine a completely 
unknown state or process, since there are exponentially 
many parameters to measure. To make tomography fea- 
sible for larger systems, we need to find a restriction that 
requires fewer measurements, ideally polynomially many. 

One way to improve the measurement, or query, com- 
plexity is to assume some prior knowledge of the process. 
For example, suppose the process was known to be one 
of a small number of unitaries, then the task is just to 
decide which. This is the approach we take here. As a 
simple example, consider being given a black box imple- 
menting an unknown Pauli matrix. By applying this to 
half a maximally entangled state, the Pauli can be identi- 
fied with one query. This is essentially superdense coding 
[7] and is explained in Section III A. Indeed, if the black 
box performed a tensor product of arbitrary Paulis on n 
qubits then it too can be identified with just one query. 
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We extend this to work for elements of the Clifford 
group (the normaliser of the Pauli group; see Definition 
1) and show that any member of the Clifford group can be 
learnt with 0(n) queries, which we show is optimal. The 
Clifford group is an important subgroup of the unitary 
group that has found uses in quantum error correction 
and fault tolerance [8-10]. 



Then generalising further, we show that elements of 
the Gottesman-Chuang hierarchy [11] (see Definition 2), 
also known as the Ck hierarchy, can also be learnt effi- 
ciently. As the level k increases, the set Ck includes more 
and more unitaries so this implies ever larger sets can be 
learnt, although the number of queries scales exponen- 
tially with k. Our methods also work if the unitary is 
known to be close to a Clifford (or any element of Ck for 
some known k) rather than exactly a Clifford. 



We also give a Clifford testing algorithm, which deter- 
mines whether an unknown unitary is close to a Clifford 
or far from every Clifford. This is an extension of the 
Pauli testing algorithm given in [12]. Indeed, our results 
are closely related to results in [12] and we use some 
of the algorithms presented there as ingredients. Our 
results can also be compared with [13], which contains 
methods to approximately learn quantum states. An- 
other related result is that of Aaronson and Gottesman 
[14], which provides a method of learning stabiliser states 
with linearly many copies. 



We only consider query complexity although, at least 
for the Clifford group results, our methods are computa- 
tionally efficient too. 



The rest of the paper is organised as follows. In Sec- 
tion II, we define the Pauli and Clifford groups and the 
Gottesman-Chuang hierarchy. In Section III we present 
our algorithm for exact learning of Clifford and Ck ele- 
ments. In Section IV we show how to find the closest 
element of Ck to an unknown unitary. In Section V we 
present our Clifford testing algorithm and then conclude 
in Section VI. 
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II. THE PAULI AND CLIFFORD GROUPS AND 
THE GOTTESMAN-CHUANG HIERARCHY 

Firstly, we define the Pauli group. The Pauli operators 
are 

ai = (o l) ° x = (l o) 

°v = ( • o*) °> = (j _°i) (!) 

which we extend to n qubits by taking all tensor products 
of these one qubit matrices. Call this set V with elements 
<jp with p e {I,x,y,z} n . We then have \V\ — 4™ and 
tr a p (jq — 2 n S pq . Also, the Pauli matrices form an orthog- 
onal basis for 2™ x 2™ matrices. Therefore we can write 
any such matrix in the Pauli basis in the form J2 P l(p)°~p- 
Then to make this into a group, the Pauli group V, we 
must include each matrix in V with phases {±1, ±i}. 
We can now define the Clifford group: 

Definition 1 (The Clifford group). The Clifford group 
is the normaliser of the Pauli group i.e. 

C = {Ue U(2 n ) : UVU^ C V}. 

Then the Gottesman-Chuang hierarchy is a generali- 
sation: 

Definition 2 (The Gottesman-Chuang hierarchy [11]). 
Let C\ be the Pauli group V . Then level Ck of the hier- 
archy is defined recursively: 

C k = {UeU(2 n ) : UVU^ C C k -!}. 

By definition, C 2 is the Clifford group C. For k > 2, 
Ck is no longer a group but contains a universal gate set, 
whereas C\ and Ci are not universal. 



III. LEARNING GOTTESMAN-CHUANG 
OPERATIONS 

Before we give our algorithm for learning Gottesman- 
Chuang operations, we present a simple method for learn- 
ing Pauli operations, which we use as the main ingredient. 



A. Learning Pauli Operations 

This is due to [12] and is in fact identical to the 
superdense-coding protocol [7]. 

Theorem 3 ([12], Proposition 20). Pauli operations can 
be identified with one query and in time 0(n). 

Proof. Apply the operator o~ p to half of the maximally en- 
tangled state IV') = 2~™/ 2 J2 { For different choices of 



a p , the resulting states are orthogonal so can be perfectly 
distinguished: 

(V| {o-p ® I) {a q ® /) |V) = 2~ n ^2(ii\a p a q ® I\jj) 
= 2- n ^2(i\a p a q \j)(i\j) 

= 2~ n J2( i M i ) 

i 

= 2~ U tl 0~pO~q 
— 3pq- 

The time complexity 0(n) comes from the preparation 
and measurement operations. □ 



B. Learning Clifford Operations 

We can now present our algorithm for learning Clif- 
ford operations to illustrate our main idea for learning 
unitarics in the Gottesman-Chuang hierarchy We will 
use the fact that knowing how a unitary acts by conju- 
gation on all elements of V identifies it uniquely (up to 
phase): 

Lemma 4. Knowing Ua p W for all a v e V uniquely 
determines U , up to global phase. 

Proof. The Pauli matrices form a basis for all 2™ x 2™ 
matrices so knowing the action of U on the Paulis is 
enough to determine the action of U on any matrix up to 
phase. The phase cannot be determined because action 
by conjugation does not reveal the phase. □ 

Now let G = {a Xi , a Zi }™ =1 where a Xi (a Zi ) is the matrix 
with a x (a z ) acting on qubit i and trivially elsewhere. 
We think of this as a set of generators for V since each 
element of V can be written as a product of elements of 
G, up to phase. Using this, knowledge of how U acts on 
elements of G is sufficient to determine the action on all 
of V: 

Lemma 5. Ua p W for any a p £ V can be calculated from 
knowledge of Ua g U' for each a g G G. 

Proof. Let a p = aa gi . . . a gm for a gi € G where a is a 
phase. Then 

VupU^ = aUa gi ...o g JJ* = aUa gi U^ . . .Ua gm Ul □ 

With these definitions and observations, we can now 
present the Clifford learning algorithm. 

Theorem 6. Given oracle access to an unknown Clifford 
operation C and its conjugate C\ C can be determined 
exactly (up to global phase) with 2n+l queries to C and 
2n to C^ . The algorithm runs in time 0(n 2 ). 
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Proof. From the definition of the Clifford group, 
Ca p C^ G V for all a p e V . Note that Ca p C^ is not nec- 
essarily a Pauli operator in V because there is a phase 
of ±1 (complex phases are not allowed because Co-pC^ 
is Hermitian). Determining which Pauli operator and 
phase for every a p would be sufficient to learn C using 
Lemma 4. But from Lemma 5, we only need to know 
Co~ g C' for each a g e G. 

Let Ca Xi = otiO- ai and Ca Zi — /3icr hi , where 
ai,(3i = ±f. Knowing just a ai and cr^ is enough to 
specify C up to a Pauli correction factor a q which gives 
the phases on and (3i. Choosing a q that anticommutes 
with a Xi flips the sign of on and similarly for a Zi . We 
now present the algorithm: 

f . Apply Ca Xi C^ and Ca Zi C^ for each i and use The- 
orem 3 to determine a ai and er^. This uses 2n 
queries to both C and CL 

2. Let C be such that C'a Xz C"< = a at and C'a z% C^ = 
a-;,; i-e. the phases are all +1. Then, choosing a 
phase for C, we can write C — C'a where 



*= n n 



(2) 



i:oti — — 1 



Then implement C'^C to determine a using The- 
orem 3. This uses one query to C. We can now 
calculate the phases ctij and 

To work out the time complexity, note that in step 1 the 
0(n) time Pauli learning algorithm is called 2n times. 
Then for step 2, the Clifford C can be implemented in 
0(n 2 ) time using for example Theorem f0.6 of [1]. □ 

We now show that this algorithm is optimal, in terms 
of number of queries, up to constant factors: 

Lemma 7. Any method of learning a Clifford gate re- 
quires at least n queries. 

Proof. Each application of the gate C can give at most 
2n bits of mutual information about C. This follows 
from the optimality of superdense coding [7]. The 
Clifford group (modulo global phase) is of size [15] 
2 „ 2 +2r l+3 rj™ =1 (4^ - 1) > 2 2 " 2 +"+ 3 . To identify an el- 
ement with to queries, we therefore need 



2 2n,n > 2 



which implies m > n. 



□ 



It is unfortunate that access to is also required, but 
we do not know a method with optimal query complexity 
that works without C< . There are however methods that 
use 0(n 2 ) queries that do not use . The result of [16] 
can be used to show that 0(n 2 ) queries to C are sufficient, 
by distinguishing the states C®I\ip) for different Cliffords 
C and where \ip) is the maximally entangled state. We 
can use Lemma 12 to show that these states are far apart 
in the distance measure used in [16], allowing us to apply 
their result. 



C. Learning Gottesman-Chuang Operations 

Theorem 6 can easily be generalised to learning any 
operation from the Ck hierarchy: 

Theorem 8. Given oracle access to an unknown op- 
eration C € Cfe and its conjugate , C can be deter- 
mined exactly (up to phase) with ( - 2 2 T ^_^ 1 queries to C 
and (2n) k - 1 to C+. 



Proof. The proof is by induction. The base case is for 
the Paulis and is proven in Theorem 3. Then, to learn 
C G Cfe+i, we assume we have a learning algorithm for 
members of Ck- Apply Co~ g C^ for each a g G G. These 
operations are elements of Ck so use the learning algo- 
rithm for Ck to determine these up to phase. Then use 
the last step of Theorem 6 to determine the phases. 

We now determine the number of queries to C and . 
Let T(k) be the number of queries to C and T'(k) the 
number of queries to . We have the recurrences 



and 



T{k+ 1) = 2nT(k) + 1 
T(l) = I 



T'(k + 1) = 2nT'(k) 
T'{2) = 2n 



(4) 



(5) 



which have solutions T(fc) = ^"J^ 1 and T'(fc) = 
{2nf- 1 (with T'(l) = 0). " □ 



IV. LEARNING UNITARIES CLOSE TO C k 
ELEMENTS 

Here we suppose that we are given a unitary that is 
known to be close to an element of Ck for some given k. 
We present a method for finding this element. But first 
we must define our distance measure. 

We would like our distance measure to not distinguish 
between unitaries that differ by just an unobservable 
global phase. We define a 'distance' D below with this 
property. However, firstly define the distance D + to be 
a normalised 2-norm distance: 



(3) Definition 9. For U\ and Ui d x d matrices, 



where \ \A\\ 2 = Vtr A^A. 

We have chosen the normalisation so that < 
D + (Ui, U2) < 1- We now define our phase invariant 'dis- 
tance': 

Definition 10. For U\ and [7 2 d x d matrices, 



D(U U U 2 ) := 



f 



\Ut®UZ -U2®U2*\\2 
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This is not a true distance since D{JJ\,U2) = does 
not imply U\ = U2, but that U\ and U2 are the same up 
to a phase so the difference is unobservable. From the 
2-norm definition, we can show: 



Lemma 11. 



D + (U U U 2 ) = 



Re 



and 



D(U 1 ,U 2 ) 



tvUiU. 



(6) 



(7) 



From this we can easily see that < D{U\,U2) < 1 
with equality if and only if U\ and U2 are orthogonal. 
Further note that by the unitary invariance of the 2- 
norm, both D and D + are unitarily invariant and from 
the triangle inequality for the 2-norm they both obey the 
triangle inequality. 

Our approximate learning method will find the unique 
closest element of Ck to U. In order to guarantee unique- 
ness, the distance must be upper bounded: 

Lemma 12. If D(U,C) < 2fc _ 1 1/2 for some C G Ck then 
C is unique up to phase. 

The proof is in Appendix A. 

Theorem 13. Given oracle access to U and W and k 
such that D(U,C) < e for some C G Ck with 



^2(1 - (2 fe - 1 e) 2 ) - 1 > 



(8) 



then C can be determined with probability at least 1 — 8 
with O ^{2n) k ~ 1 log (2 " + ] )fc ~ 1 ) queries. 

Proof. By Lemma 12, C is unique up to phase. We now 
prove the Theorem by induction. 

For k = 1, use Proposition 21 of [12] to learn the closest 
Pauli operator. This works by repeating the Pauli learn- 
ing method Theorem 3 and taking the majority vote. 
This uses O (-^ log |) queries to succeed with probabil- 
ity at least 1—8. 

Now for the inductive step. Assume we have a learn- 
ing algorithm for level k. Then for C G Cfe+i, let 
Ccr Si C^ = C 9i for o gi G G. By Lemma 16, we have 
D(Ua gi W , C 9i ) < 2e. Use the learning algorithm for 
level k to determine C 9i up to phase for all i. Then to 
find the phases we use the same method as before: im- 
plement any C with C'o- gi C^ = ±C 9i for any (known) 
choice of phase. Then C' = Ca q for some Pauli opera- 
tor a q . We can determine a q by implementing C'^U and 
using the k — 1 learning algorithm since 



D(C^U,a q ) = D{U,Ca q ) 
= D(U,C) < e. 



Now we calculate the success probabilities and number 
of queries. There are 2n + 1 calls to the algorithm at 
lower levels, which all succeed with probability at least 
1 — 8. So at this level the success probability is at least 
1 — (2n + 1)8. So to succeed with probability at least 
1 — 8 we must replace 8 with 8/(2n+ 1). Then the overall 
number of queries is 



2n- 0(>n)- log ^ 



= 0^(2n) fc log 



(2n + l) k 



(10) 



□ 



(9) 



We remark that there is only O(fclogn) overhead (for 
constant e' and 8) over the exact learning algorithm of 
Theorem 8. 



V. CLIFFORD TESTING 

Here we present an efficient algorithm to determine 
whether an unknown unitary operation is close to a Clif- 
ford or far from every Clifford. Whereas the previous 
results allow us to find the Clifford operator close to the 
given black box unitary, in this section we are concerned 
with determining how far the given unitary is from any 
Clifford. We do not measure this directly, but provide 
an algorithm of low query complexity that decides if the 
given unitary is close to a Clifford or far from all. This 
type of algorithm is known in computer science as a prop- 
erty testing algorithm and has many applications, includ- 
ing the theory of probabilistically checkable proofs [17]. 
The result in this section could be extended to work for 
any level of the Gottesman-Chuang hierarchy although 
for simplicity we only present the version for Cliffords. 

The key ingredient to our method will be a way of 
estimating the Pauli coefficients: 

Lemma 14 (Lemma 23 of [12]). For any p G {/, x, y, z} n 

and unitary U, \~f{p)\ = ■^\trUa p \ can be estimated 

to within ±77 with probability 1 — 8 using O log 
queries. 

This is a generalisation of Theorem 3 and the method 
is similar. Instead of there being only one possible out- 
come, now the probability of obtaining the outcome cor- 
responding to a p is estimated. This probability is equal 

to | 7 (p)| 2 . 

Theorem 15. Given oracle access to U and W with the 
promise that for < e < 1 either 

a) CIOSE: there exists C G C such that D(U,C) < 
—4 — or 

V32n 

b) FAR: for all C G C, D(U,C) > e and there exists 
C G C such that D(U, C) < 1/3 
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holds then there is a O ^frlog algorithm that deter- 
mines which with probability at least 1 — S. 

Proof. In both cases, we have that D(U,C) < 1/3 for 
some C, which ensures that C is unique (using Lemma 
12, since | < and can be found using Theorem 13 

with O (nlog j) queries. Then the algorithm is: 

1. For each a g € G, measure the Pauli co- 
efficient of Co-gC^ in Uo-gW (i.e. measure 
|tr Ua g WCa g C^\ /2") to precision using 
Lemma 14. 

2. If all the coefficients are found to have modulus at 
least 1— then output CLOSE else output FAR. 

This works because, for the two possibilities CLOSE and 
FAR: 



a) Using Lemma 16, D(U, C) < 
all a p e V, 

D{Ua p U\Ca p C r ) < 



/32n 



implies that for 



2e 



32n 



(11) 



Since we will only apply Ua g U' f for a g G G we 
restrict this to only the generators to find that for 
all G g e G, 



D(Ua g U\C(7 g C ] ) < 



2e 



giving 



tr Ua g WCa g C^ 



2™ 



> 1 



32n 



e 



(12) 



(13) 



for every generator a g . We need a bound on the 
non-squared coefficients, which follows directly: 

r 2 



tr U a g WC a g C 



>!--%. (14) 
~ 8n 2 V ' 



Therefore when measuring the coefficients to preci- 

2 a 2 

sion jg^2, all results will give at least 1 — y§^2- 

b) Using the contrapositive of Lemma 17, D(U, C) > e 
implies that there exists a p <E V such that 



D + {Ua p U\C(j p C^) > e. 



(15) 



Using the contrapositive of Lemma 18 this in turn 
implies there exists a g e G such that 



D+(U*gU\C<JgCl)>-, 



(16) 



which means that for at least one a g € G, Ua g U* 
will have a small overlap with Co- g C^ i.e. there ex- 
ists a g € G such that 



tr Uo-gWCo-gC^ 



The C returned by the application of Theorem 13 
is such that trUa g U^Ca g C^ is positive, which jus- 
tifies inserting the absolute value signs above when 
using D + rather than D. This implies that at least 
one coefficient will be found to be less than 1 — 

2 

when measuring to precision □ 



VI. CONCLUSIONS AND FURTHER WORK 

We have shown how to exactly identify an unknown 
Clifford operator in O(n) queries, which we show is op- 
timal. This is then extended to cover elements of the 
Ck hierarchy and for unitaries that are only known to be 
close to Ck operations. The key to the Clifford learning 
algorithm is to apply Co- p C^ and then find the resulting 
Pauli operator. 

A way of extending this idea could be to learn uni- 
taries from larger sets. Suppose S is a set of unitaries 
with the property that for every S € 5, So-pS^ is a linear 
combination of a constant number of Paulis. Then S can 
be learnt in the same way as above, using the quantum 
Goldreich-Levin algorithm of [12], which can efficiently 
find which Paulis have large overlap with an input uni- 
tary. However, we have not been able to find interesting 
sets S other than the Clifford group with this property. 

We also presented a Clifford testing algorithm, which 
determines whether a given black-box unitary is close to 
a Clifford or far from every Clifford. This can be seen 
as a quantum generalisation of quadratic testing, just as 
Pauli testing can be seen as a quantum generalisation of 
linearity testing. Property testing of this form is used 
to prove the PCP theorem [17] so these quantum testing 
results could potentially be useful in proving a quantum 
PCP theorem. It would also be interesting to strengthen 
the testing method in Theorem 15 to remove the 0(l/n) 
difference between the close and far conditions. 

Finally, it would be interesting to sec if it is possible to 
remove the requirement to have access to W . However, 
using both U and is the key to our method so we do 
not know if a method without is possible with low 
query complexity. 
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APPENDIX A: PROOF OF LEMMA 12 

Proof of Lemma 12. The proof is by induction. The base 
case is for k = 1 when we have the Pauli group. Without 
loss of generality, assume C is a Pauli operator with no 
phase. Let C — a p . 

Expand U in the Pauli basis: 



(Al) 



Since U is unitary, we have J2 q It | 2 = 1- By Lemma 
11, 



D(U, a v f = 1 - 



tr o~ p U 



2 n 



which implies 



|7(P)| 2 > 1 



(A2) 



(A3) 



Now, suppose for contradiction that there exists o Pl ^ 
a P2 with D{U,o pi ) < e and D(U,a p2 ) < e. Then by the 
above, |7(pi)| 2 , | T (P2 ) 1 2 > 1 — e 2 - But there is also the 
constraint |7(pi)| 2 + |t(P2) | 2 < 1 which combined give 



e > 



V2 



(A4) 



which is false by assumption. This implies a Pl — cr P2 , 
which proves the base case. 

To prove the inductive step, again assume for contra- 
diction that there exist C\,C 2 € Ck+i with C\ ^ C 2 
and D(U,Ci) < e and D(U,C 2 ) < e. Then there exists 



cT g G G with 



C\G g C\ =: C\ g 7^ C 2g := C 2 a g C\. 



(A5) 



Here, C lg , C 2g € C k . 

Using Lemma 16, D{Ua g U\Ci g ) < 2e and 
D(U<j g U\C 2 g) < 2e. 

Now there are two cases. Firstly, suppose we can 
choose a g such that C\o~ g C\ ^ ±C 2 o~ g C 2 . Then C\ g and 
C 2g are not equivalent up to phase so, using the inductive 
hypothesis, we must have 



or 



2e > 



e > 



2 fc-l/2 



1 



2 (fe+l)-l/2 



(A6) 



(A7) 



which is again false by assumption. 

For the other case, C\a g C\ = ±C 2 o g C\ for all a g € G. 
This implies that C 2 — C\o q for some Pauli a q ^ /. Then 
we have 



£>(C/,d)<e 
D(U,CKT g ) < e 



(A8) 



which by unitary invariance gives 

D(C\U,I) < e 
D(C\U,a q )<e. 



(A9) 



But we proved that this is impossible in this range of e 
in the k = 1 proof above. □ 



APPENDIX B: MISCELLANEOUS LEMMAS 

Here we prove some miscellaneous lemmas used earlier 
in the paper. 

The first lemma says that for two close operators U\ 
and U 2 , U\u p U\ is close to U 2 o p \j\ for all Paulis a p : 

Lemma 16. If D(U\, U 2 ) < 8 then for all a p e V, 

D{Ukt p u1,U 2 <t p U%) < 26. 

Proof. Let U\ = VU 2 and U 2p — U 2 a p \j\. Then we sim- 
ply apply the triangle inequality for D and unitary in- 
variance: 

D^VpUl U 2 o p u\) = D{VU 2p V\ U 2p ) 
= D(VU 2p ,U 2p V) 
< D(VU 2p , U 2p ) + D(U 2p , U 2p V) 
= D(V,I)+D(I,V) 
= 2D(U U U 2 ). □ 

The next lemma is a converse to this: 
Lemma 17. If for all a p e V 

D + (U 1( T p Ul U 2 u p U\) < 5 (Bl) 

then 

D(U U U 2 )<5. (B2) 

Proof. If D+{Ui<T p Ul,U 2 o- p Ul) < 5 then 
TprRctYUio-pUlU^pUt, > 1 - S 2 . Since this is 
true for all a pj we can take the average of this over the 
whole of V and use the fact that for any d x d matrix A 
W^a- ^f a v-^ a p — 2^" tr ^ (the Paulis are a 1-design) 
to find 

Retr U, tr u\U 2 ^j U\>1- S 2 (B3) 



which simplified gives 



trLW] 



> 1-S 2 



giving the desired result. 



(B4) 



□ 



Now we show how to go from distances for just the 
generators G to distances for the whole of V: 
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Lemma 18. If for all a g e G Proof. The proof is by induction on the number of gener- 

ators required to make a pi using the triangle inequality 
D+(U 1 a g Ul,U 2 a g U^) < S (B5) for £>+. □ 

then for all a p € V 

D+(U l0 - p Ul, U 2 (J p Ul) < 2nS (B6) 
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